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Many research groups are starting to pay serious attention to the problem of 
fractional-order circuits. In this paper, a new approach to designing fractional 
order low-pass filter (FOLPF) is presented. Finding a rational approximation of 
the fractional Laplace operator s“ is a crucial step in the design of fractional 
order filters. A comparative study of the most widely used approximation tech- 
niques named continued fraction expansion (CFE) method and Biquadratic Ap- 
proximation (RE) method is performed. Then the transfer function of the pro- 
posed FOLPF is calculated. Using operational amplifier, the proposed filter is 
synthesized. The proposed circuit is simulated using Texas instruments TINA 
software. The results obtained outperform the existing methods. 
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1. INTRODUCTION 

The simulation of the physical environment can make use of an extension of classic integral and dif- 
ferential equations known as fractional order differential and integral equations |I]. The fractional differential 
equations could not be solved using any of the methods that are available, so the integer-order models were used 
instead. Several different approaches are now being utilised in the process of developing applications predi- 
cated on fractional derivatives and integration. It is possible to improve the characterization of a wide variety 
of dynamic systems by employing non-integer order dynamic models that are founded on fractional calculus, 
differentiation, or integration. Calculus in its classical form is built on the foundation of differentiation and 
integration of integer orders [I]. Fractional order systems makes use of fractional order differentiation and 
integration as constituent elements. Recent applications of fractional order systems in electrical engineering is 
an appreciable thing. A fractional order element (also called as fractance device) is designated with s“. 

The design of fractional order filters requires a number of steps, one of which is the realization of 
the fractance device. Krishna and Reddy discussed a realization of a fractance device that is based on 
a continuing fraction growth. The realization of a fractional impedance of order 1/2 is made possible by 
Hilbert’s structure, which is also put to use in the construction of fractional order controllers [3]. Many of 
the electrical systems make use of something called a fractional order filter. The literature contains many 
examples of very well built low pass, high pass, band pass, band reject, and all pass fractional order filters, 
as well as first and second order pass filters. They are built to operate in both the current mode and the 
voltage mode simultaneously. In the design of a fractional-order filter, some of the many active elements 
that are commercially available and can be used include operational amplifiers, differential difference current 
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conveyors (DDCCs), current feedback operational amplifiers (CFOAs), second generation current conveyors 
(CCIls), operational transconductance amplifiers (OTAs), and current amplifiers (CAs) [4]. To generate a 
transfer function that is analogous to that of a fractional-order filter, one may employ any one of a number 
of different topologies. The topologies follow the leader feedback (FLF), inverse follow the leader feedback 
(IFLF), or single amp biquad (SAB) in conjunction with a first-order filter are the ones that are utilised the 
most frequently. Even if there are a variety of topologies available, the end user must make use of a specific 
topology in order to satisfy the design requirements. 

Doréak et al. discussed the design and electronic implementation of the fractional-order controller 
and controlled system. This system is based on the Laplace transfer function of the corresponding electronic 
circuit being equivalent to the Laplace transfer function of the original fractional order systems. Mahmood et al. 
[6] purposed to synthesise a fractance circuit in an electronic circuit to achieve fractional order derivatives. Two 
different types of fractances were synthesised and linked in an operational amplifier as passive RC circuits to 
achieve the realization. Iqbal and Shekh [/7] have studied about different types of approximations of fractional 
order systems. Realization of fractional order filters at microwave region is studied in [8]. Tewari and Arya [9] 
have designed and performed the PSPICE based simulation of operational transresistance amplifier (OTRA) 
based fractional low pass filter of order 1.5. Resistor less realization of generalized filters of fractional order 
using OTAs as active elements is realized in [10]. Operational amplifier based design of fractional order 
controllers is studied using simulation tools [IT]. Sacu and Alci used Cadence-PSpice to build and simulate 
simple OTA-based fractional order low-pass and high-pass filters of order (n + a), where 0 < a < 1 and 
n > 1. For any order 0 < a < 1, the fractional-order low-pass filter based on a single fractance element 
is represented in by using optimal integer-order transfer function approximations. Initially, integer-order 
filter coefficients for the FLF are found directly in the range of 0.01 to 0.99 with a step size of 0.01. This is done 
by employing the colliding bodies optimisation technique. A novel rational approximation for the fractional 
order operator is proposed in [14]. Signal flow graphs synthesise OTA based fractional element in [15]. Theory, 
simulation, and approximation function match. Genetic algorithm is used in the synthesis and optimization of 
fractional order elements [16]. Different passive realization of fractional order capacitors studied in [17]. 

In order to design the fractional order low-pass filter (FOLPF), it is necessary to approximate the 
operator to a certain degree on a numerical scale. These methodologies for approximation were proposed by 
the researchers at various points in time throughout the study. Based on a general filter topology, the work that 
was done in introduces fractional order low-pass, high-pass, band-pass, and band-reject filter responses at 
the same circuit utilising low-voltage active elements. Simulation of the designed circuits is performed with 
Spice software and makes use of the properties of 0.35 um Taiwan semiconductor manufacturing company 
(TSMC) complementary metal-oxide semiconductor (CMOS) technology. This article presents the results of 
a comparative examination of the many rational approximation approaches that are currently available. Yüce 
and Tan proposed a basic way for generating fractional-order filters with transfer functions by making use 
of Laplace operators of a variety of fractional orders. Mishra et al. have come up with a proposal for a 
fractional order butterworth lowpass filter that makes use of differential voltage current conveyor. 

A discussion of a proven example of the realization of a fractional order element through the use of the 
partial fraction expansion method in [2I], [22]. Theirs work compares two different approaches to design by 
imitating skin tissue and using AgCI electrodes in order to test and evaluate electrical bio-impedance circuits 
and systems [23]. The models are built on fractional-order elements, implemented with active components, 
and may reflect bio-impedance characteristics up to 10 kHz. A fractional-order passive RC low-pass filter 
is investigated in [24]. The progression of behavior in the temporal realm was revealed through fractional 
orders. MATLAB’s output consisted of low-order fractional-order filters. The innovative method of curve 
fitting provides a close approximation of this function. Using a universal fractional-order transfer function, 
this method provides an explanation for all of the many types of fractional-order filters. A fractional-order 
element analogue modelling system is illustrated in [25], which may be found here. The traditional circuit 
theory is applied here. Mijat et al. discusses a number of different processes that can be used to make 
the fractional element. A fractional change occurs in the simulated frequency response. A single OTRA is 
utilised in the construction of an active current mode fractional order analogue filter. The 1.5-order filter that 
has been proposed. The findings of the PSpice simulation matched the theoretical predictions. The RC-RC 
decomposition was utilised in the filter that was manufactured by OTRA. Prommee et al. demonstrate 
the application of fractional order devices in the processing of bio-medical signals employing OTAs. These 
amplifiers work in conjunction with the fractional order devices. Research by Krishna provides an in- 
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depth analysis of the many developments that have been made in the realm of realising the fractance device. 
A variety of different comparable circuit representations of fractional order elements are given by [29], [80]. 
Narayan et al. present a realization of a fractional order bandpass filter by employing a reconfigurable 
device. The paper is organised as follows: different rational approximation methods are discussed in section 2, 
section 3 deals with the FOLPF, section 4 discusses the implementation of the filter, the finalized circuit and 
the corresponding results are presented in section 5, and finally conclusions are drawn in section 6. 


2. RATIONAL APPROXIMATION METHODS 

Given that the fractional order system is in fact an infinite order system, it is preferable to approximate 
lengths with finite values. Both the Laplace and the z domains can be used to make an approximation of the 
fractional order operator known as s“. Approximations that fall within the s-domain umbrella include the 
oustaloup approach, the mastudas method, and the least square method. Each approach has a number of 
advantages and disadvantages. The continued fraction expansion (CFE) method and Reyad Elkhazali (RE) 
approximation are both taken into consideration in this article as potential methods of approximation. Because 
of the complexity of the hardware, only approximation up to the second order is taken into consideration. It 
can be written as (1) [4], which is the second order approximation of s®. 


2 
a _ aos +aıs+az 
a= a28? +a1ı s+a0 (1) 


where ao, @1, az are the filter coefficients. 


2.1. Continued fraction expansion method 
The CFE of (1 + x)® is used to characterise this technique [4], [32]. 


(1 $ x) — 1 ag (1+a)s (1—a)xr (2+a)x (2—a)x (2) 


I= lF 24 3+ 2+ 5+... 


when x = s — 1 is substituted, the above CFE converges down the negative real axis from s = —oo tos = 0 
in the finite complex s-plane. Considering the first four terms of the expansion a second order approximation 
will be derived with the coefficients as (3): 


ag = a? + 3Ba+2 
a, = 8— 2d? (3) 
a2 =a? — 3a +2 


2.2. Biquadratic Approximation method 

Khazali et al. presented this approximation in 2019. This method is designated as RE to designate 
the name of the scientist. A biquadratic approximation approach can be used to estimate a fractance operator. 
The rational approximation will take the shape of a cascade connection of biquadratic transfer functions as (4): 


(G) -ia (2) - i z o 


i=1 Di( aay) 


where 


wi = Center Frequency 


wg = Geometric Mean (5) 
The remaining center frequencies may be computed using a recursive formula using w; as in (6): 
Ww = an a: i= 2,3,4,...,n. (6) 
Wg is calculated by solving in (7): 
agagny* + a1 (az — a9)7? + (at — a3 — a)n? + a1 (a2 — ao)Y + aoaan = 0 (7) 


where 1) = tan(%F ). The expression for the quadratic will be (8): 
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The values of the coefficients ap, a1, a2 are calculated in (9): 
ao =a% +2a+1 
ag =Q% — 2a +1 
aı = (az — ao) tan (C=) 


After calculating the coefficients, the rational approximation for s~ may be determined. The second order 
rational approximations obtained for different values of a is as shown in Table 1. 


3. 


Magnitude[dB] 


Table 1. Rational approximation of second order 


a Hre(s) Hor x(s) 
0.1 1.99487 +5.0828s+1.594 2.31s7+7.98s+1.71 
g 1.594s2+5.082s+1.994 I.71s2+7.98s+2.31 
0.2 2.1258?+5.051s+1.325 2.6457+7.92s+1.44 
g 1.32582+5.051s+2.125 1.44s2+7.92s+2.64 
03 2.29787+4.998s+1.097 2.99s?+7.82s+1.19 
: 1.097s2+4.998s+2.297 1.19s2-+7.828+2.99 
0.4 2493s? +4.924s+0.8931 3.3657+7.68s+0.96 
i 0.893152 +4.924s+2.493 0.96s2+7.685+3.36 
0.5 2.707s?+4.828s+0.7071 3.75s? +7.5s+0.75 
s 0.707182 -+4.828s+2.707 0.7582 +7.88+3.75 
0.6 2.93652 +4.71s+0.536 4.16s7+7.288+0.56 
i 0.53652 +4.71s+2.936 0.56s2+7.28s+4.16 
07 3.17987 +4.569s+0.3791 4.59s?+7.02s+0.39 
j 0.379152 +4.569s+3.179 0.39s2+7.02s+4.59 
0.8 3.4375? +4.404s+0.2365 5.04s?+6.72s+0.24 
: 0.236582 +4.404s+3.347 0.24s246.72s+5.04 
0.9 3.71s7+4.215s+0.1095 5.51s7+6.38s+0.11 
d 0.109552 +4.215s+3.71 0.11s2+6.38s+5.51 


A comparison of the two rational approximations (for œ = 0.4) is as shown in Figures | and 2. The 
phase response offered by the RE approximation is good as compared to CFE based method. From the magni- 
tude response it is observed that the CFE based method produces good results for larger range of frequencies 
as compared to RE method. Keeping in view of the better linearity in phase angle, RE method is considered 
for the implementation. 
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Figure 1. Comparison for magnitude response of 
different approximation methods 


FRACTIONAL ORDER LOW-PASS FILTER 


Figure 2. Comparison for phase response of 
different approximation methods 


Fractional order filters can meet the exact design criteria, but traditional integer-order filters were 
unable to do so. Comparing the design parameters to the integer-order filters, the fractional notion increases 
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the level of safety a little amount. Fractional order filters make it simple to meet the requirements for the circuit 
design and tuning parameters. Only fractional order filters make it feasible to vary and modify parameters like 
the roll-off frequency to any desired slope. Thus, compared to integer-order filters, fractional order filters offer 
greater design freedom. 

A standard filter and a fractional-order filter vary in that additional parameters are included in the term 
that determines the slope of attenuation of a given transfer function. Therefore, the equation for a fractional- 
order filter’s attenuation of transfer function is 20(n + a) (dB/decade). Where n is an unsigned integer number, 
often between 1 and 10, and a is defined as a real number in the range [0 1]. The transfer function of FOLPF 
is given by (10) [24]: 


LP _ k 
Anta = S87 tha)t ks (10) 
where n is an integer and 0 < a < 1. The following are the values of k’s as (11): 
ky =1 


k2 = 1.17960? + 0.1670 + 0.21735 (11) 
k3 = 0.19295a + 0.81369 


Substituting the second order approximation of the fractional order operator in the transfer function of the 
low-pass filter, the equation becomes as (12): 


HLP = kı 


ai a s2 asta 
MEES a e ) (on tha) ths (12) 


After rearranging the terms, the equation becomes as (13): 


ky (ays? +aıs + ao) 


HEP = 13 
n+a ags” t? 4 azs” t! L ags” ae s2 (k3a2 te kao) ( ) 
+s (k3a1 + kga1) = (k3ao + k2a2) 
For n = 1, in (13) simplifies to (14): 
EP kı (azs?-+a1s+a0) 
Avie ~~ ao(s8+c1 s2+c28+c3) (14) 


where, 


cı = (k3a2 + k2ao + a1) /ao 
C2 = (kza1 -= kya, "E a2) / a0 (15) 
C3 = (k3ao T k2a2)/ao 


The equations indicate that the values of the filter coefficients are dependent upon the value of frac- 
tional order œ. Table 2 depict the values of filter coefficients for different value of fractional order. The transfer 
functions of the FOLPF using different rational approximation techniques (for different values of œ) as tabu- 
lated in Table 3. 


Table 2. Values of filter co-efficients for different values of a 

Qa ao ay, a2 Cl C2 C3 

0.1 2.31 7.98 1.71 3.497325 0.9053559 0.04038944 
0.2 264 7.92 1.44 3.101633 0.9029765 0.08254491 
0.3 2.99 7.82 1.19 2.794686  0.9580753 0.1200770 
0.4 3.36 7.68 0.96 2.563302 1.046208 0.1501903 
0.5 3.75 7.5 0.75 2.397695 1.149750 0.1721550 
0.6 416 7.28 0.56 2.29044 1.255711 0.1864237 
0.7 4.59 702 0.39 2.235792 1.354332 0.1941091 
0.8 5.04 6.72 0.24 2.229228 1.438158 0.1966716 
0.9 5.51 6.38 0.11 2.267138 1.501410 0.1957304 
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Table 3. Transfer functions of FOLPF 


a Hrpore(s) HrpRE(s) 
0.1 1.71s2+7.98s+2.31 1.5948” +5.0820s+1.9940 

i 2.31s3+8.0788s?+2.0914s+0.0933 1.9940s3+5.1696s? +1.8369s-+0.0839 
0.2 1.44s? +7.925s+2.64 1.325s?+5.051s+2.1250 

i 2.6453 +8.1883s2+2.3839s+0.2179 2.125053 +5.2734s2+1.9269s+0.1888 
03 1.19s?+7.82s+2.99 1.0975? +4.998s+2.297 

f 2.99s3+8.3561s2+2.8646s+0.359 2.29753 +5.4204s2 +2.1973s-+0.3044 
0.4 0.96s?+7.68s+3.36 8.931s7+4.9240s+2.493 

: 3.3689 +8.612752+3.5153s+0.5046 2.49383 +5.63s82+2.5314s+0.4206 
05 0.7587 +7.58s+3.75 7.07187 +4.8280s+2.707 

i 3.75s3+8.9914s2+4.3116s+0.6456 2.70753 +5.9205s7+2.9998s+0.5287 
0.6 0.56s7+7.285+4.16 0.53687 +4.71s+2.9360 

` 4.16s3+9.52825s2+5.2238s+0.7755 2.936053 +6.3130s2+3.5533s+0.6212 
0.7 0.39s7+7.02s+4.59 3.791s?+4.569s+3.179 

: 4.5985 +10.2623s82+6.2164s+0.8910 3.17983 +6.8293s82+4.1712s+0.6928 
0.8 0.24s?+6.72s+5.04 0.23655? +4.4040s+3.437 

i 5.0483 +11.2353s2+7.2483s-+0.9912 3.43783 +7.4944s2+4.8294s+0.7407 
0.9 0.11s?+6.38s+5.51 0.10955? +4.215s+3.71 

i 5.51s3412.4919s248.2728s+1.0785 3.7183+8.3364s2+5.5023s+0.7653 


4. IMPLEMENTATION OF THE PROPOSED FILTER 


In order to implement the transfer functions derived in the previous sections, FLF topology is used. 
In this paper, feed-forward path based FLF structure is used for implementation. The outputs are available in 
terms of s” Vout, STI Vouts «+++ , SVout, Vout- These outputs are used in a feed-forward form to get nt” order 
filter with arbitrarily transmission zeros as shown in Figure 3. In this structure, the outputs from the lossless 
integrators are multiplied by coefficients and then added together to get final output [33], [34]. 


Vi. _ 1 
Vag T Hils) ae pka kp 1e Fkh (16) 
Tys = Ho(s) = aos” + as”! + ags”—?.....dn_18 + an (17) 


‘out = H(s) = Hy(s)Ho(s) 
— aos” pars” Ip ags P.n. Qn—18+an (18) 
s @+kys"—-l+kos"—-2..... Kkn—18stkn 


The above transfer function can be realized using only lossless integrators and multipliers. 


Vout: Vis — — _ As?+Bs+C 
Vin Vi Vin Hı H = s3+Ds2?+Es+F (19) 


Now splitting the transfer function. 
Hı = Yet = As? + Bs +C 


=> Vout + (—As?V, — BsV, — CV) = 0 (20) 
> Vout +r ( s? Vi WA sVı T/B Vi ve) = 0 


_ MN _ 1 
H2 = yh = sypstEeiF (21) 


The above equation can be rearranged as (22): 


Vin — (s° + Es)V, — (Ds? + F)V, =0 


= Vin~ (+ s1) Vi (217o + afr) Vi =0 a 


The operational amplifier based realization valid for any fractional order a is as shown in Figure 4. 
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Figure 4. Synthesized circuit for any value of a 


5. RESULTS 
To illustrate the methodology, a value of a = 0.7 is considered. In (23) is the transfer function of RE 
approximation is taken into consideration. 


LP _ 3.791s744.569s+3.179 
Hi7 (s) = 3.179s3+6.8293s2+4.1712s+0.6928 (23) 


The circuit realized is as shown in Figure 5. If frequency denormalization by a factor of 20k rad/sec and 
impedance scaling by a value of 10000, then all resistances and capacitance values will be practicable. The 
same circuit is used for the simulation purpose using TINA software. For the simulation purpose Texas instru- 
ments TINA software is used. The corresponding results of transient response and bode plot are as shown in 
Figures 6 and 7 respectively. 
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Figure 6. Transient response 
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6. CONCLUSION 

This paper presents a novel realization of FOLPF using different approximation techniques. Initially, 
a comparative study of the existing approximation techniques for fractional order device is carried out. Then 
synthesis of the FOLPF is carried out using FLF topology. The proposed method is straight forward and makes 
use of an operational amplifier for the realization. This method makes use of only lossless integrators and 
amplifiers. A 100 mHz and 1 Vpp sinusoidal signal is chosen as input. LM301A operational amplifier is 
chosen for simulation purpose. It is observed that there is a phase shift of 180 between input and output. The 
linear range of operation can be up to 100 uHz to 2 KHz. So, the procedure suggested can be used for the 
realization of fractional order filters. 
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